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^ \ Abstract 
00 , 

, We use the available data on e+e hadrons to determine the value of the QED 

' coupling at the Z pole. The inclusion of recent preliminary BES-II data in the analysis 

is found to reduce the ambiguity in the determination of the coupling, particularly that 
arising from the use of either exclusive or inclusive data in the energy region -^s ^ 2 GeV, 
- ' in favour of a larger value of a{M'^)~^ = 128.978 it 0.027. As a consequence the predicted 

value of the mass of the (Standard Model) Higgs boson is increased, so that the preferred 



^ , value is close to the LEP2 bound. 



X 



^ [ The value of the QED coupling at the Z pole, a(M|), is the poorest known of the three 

parameters (G^?, M^, a(M|)) which define the standard electroweak model. Indeed it is the 
precision to which we know a(M|) which limits the accuracy of the indirect prediction of the 
mass Mh of the (Standard Model) Higgs boson |[l|, 0. In fact the predicted allowed domain of 
Mh appears to be not too far from being in conflict with the direct LEP2 bound on the mass. 
However the quoted mass range for the Higgs does not allow for uncertainty in the value of 
a(M|) itself, and so the conflict is less severe than implied by the profiles versus Mh |1|, 0- 
Clearly a more precise determination of a(M|) is especially important. 

The value of a(M|) is obtained from 

a-^ = a{0)-^ = 137.03599976(50) (1) 



using the relation 



ais) 



(l - A«,ep(5) - Aaa(.) - Aa*°P(.)) a~\ (2) 
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where the leptonic contribution to the running of the a is known to 3 loops 



Aaiep(M|) = 314.98 x 10-^ (3) 

From now on we omit the superscript on Aahad and assume that it corresponds to five 
flavours. We will include the contribution of the sixth flavour, Aa;*°P(M|) = —0.76 x 10^^, at 
the end. To determine the hadronic contribution we need to evaluate 

/ N "-5 /""^ R(s')ds' 

at s = M|, where R = cr(e+e^ hadrons)/cr(e"'"e~ 

The early determinations of Aa^adiM^) made maximum use of the e~^e~ measurements 
of R{s), using the sum of the exclusive hadronic channels (e"^e^ —>■ 27t,37t, . . . KK, . . .) for 
a/s < 1.5 GeV and the inclusive measurement of a{e^e^ hadrons) at larger energies, culmi- 
nating in the analysis of Eidelman and Jegerlehner |^ which gave 

Aahad(M|) = (280.4 ±6.4) x 10"^ (5) 

Other similar determinations, which gave compatible results, can be found in Refs. [§, |§, 0]) 
where the latter analysis incorporated information from r (hadrons +i^r) decays. However, 
in the last few years the determinations of Aa^ad from (^ have relied more and more on 
theoretical input. First, perturbative QCD was used to better describe R{s) (for ^/s > 3 GeV) 
in continuum energy regions above the resonances up to the next flavour threshold Then, 
encouraged by the success of perturbative QCD to describe r decay, it was used down to 
1.8 GeV, across a region with sparse data on R{s), giving 

Aahad(M|) = (277.5 ± 1.7) x 10-^ (6) 

where ±10~^ comes from the uncertainty in perturbative QCD. A recent analysis |10], which 
follows the more conservative perturbative QCD input of ||^, flnds 

Aahad(M|) = (274.18 ±2.52) x 10"^ (inclusive) (7) 
Aahad(M|) = (276.97 ±2.90) x 10"^ (exclusive) (8) 

according to whether R is evaluated in the interval 1.46 < ^/s < 2.125 GeV using either the 
inclusive data for e~^e~ hadrons or the sum of the data for the exclusive channels. The 
difference in the input in this region can be seen in Fig. 1. 

Other recent calculations of a(M|), based on (^, have used the analytic behaviour in the 
complex s plane in attempts to reduce the error coming from the measurements of R, see, for 



example, Refs. |]TT|-[jT3[. These methods have been critically reviewed by Jegerlehner [Q. In 
this work he showed that a space-like evaluation of (^ at s = —(2.5 GeV)^, followed by an 



analytic continuation |]T5[ to s = — M|, and then around the semicircle to s = M|, gives 

Aahad(M|) = (277.82 ± 2.54) x 10-^ (9) 

in agreement with the direct determination of (^ at s = M|. He concludes that no appreciable 
reduction of the error, arising from the data, is obtained by the analytic methods. The reduction 



2 



in the error in in comparison to that in is due to the use of perturbative QCD to 

evaluate R in energy regions where it is expected to be rehable. 

We see that the values (|) and (|^) are in good agreement with each other, but that the 
determination (|^ using inclusive data below 2.1 GeV is significantly lower. This discrepancy 
has existed for some time, although most recent determinations have relied on the sum of the 



exclusive data. However BES-II 1999 data have recently been presented [jT6[ which, although 
preliminary, already clarify the behaviour of R{s) in the critical 2 < ^/s < 3 GeV domain. 
These data update the previous, published, BES-II measurements [jl^. They are shown by the 
solid circular data points in Fig. 1. First, we see that in the region ^/s ~ 2.5 GeV the new mea- 
surements are closer to the expectations of perturbative QCD than the previous measurements. 
Second, at ^/s = 2 — 2.2 GeV they join reasonably smoothly with the inclusive measurements 
of the 772 collaboration [|18|, whereas they lie in the region of, or below, the lower limit of the 
determination of R from the sum of the exclusive channels. 

Here we re-evaluate @ at s = M|, incorporating these new BES-II data. To be specific 
we use the dashed curve for ^/s > 1.46 GeV for the so-called "inclusive" determination of 
Aahad(Af|), and the average of the bounds for ^/s < 1.9 GeV for the "exclusive" determination. 
In addition we use the data, and follow the method, described in p|, |10|- The errors on the input 



'data' values of R{s') are calculated using a correlated minimization to combine different 
data sets, as described in detail in Ref. 0. 

The data from ^/s = 2.8 GeV up to the DD charm meson threshold are in good agreement^] 
with perturbative QCD, apart, of course, in J/ip and tp' resonance regions. In the energy 
intervals where perturbative QCD is believed to be valid (2.8 < ^/s < 3.74 GeV and y/s > 
5 GeV) we use both the two-loop expression with the quark mass explicity included and the 
massless three-loop expression calculated in the MS renormalization scheme^. In addition 
to the above, we include the perturbative QCD errorQ coming from varying mc, mb, Mz within 
the uncertainties quoted in the p2|, as^M^) = 0.119 ± 0.002 and varying the scale as{cs) in 
the range 0.25 < c < 4. The contributions to Aahadl^l) (il) froiii the various s' intervals 
are listed in Table 1. 

The recent Novosibirsk measurements have considerably improved the exclusive data below 
1.4 GeV, see, for example, Refs. For instance, the measurements of e~^e^ —>■ tt+tt" 

p3| are of much improved precision such that, when taken together with chiral perturbation 
theory estimates, we find that the error on the 27r contribution to Acthad for ^/s < 0.96 GeV 
is reduced to about ±0.5 x 10~^. We checked that our results for the contributions from the 
exclusive channels were in agreement with the detailed table of results given in Ref. , if we 
were to omit the new Novosibirsk data. 

From Table 1 we see that the largest uncertainties in Aa^ad occur in the 1.46 < < 
2.8 GeV interval, although the introduction of the new BES-II data |jl6| have led to a significant 



improvement. It is interesting to note that if we were to assume that perturbative QCD gave 



^The MARK I data |19 are normalized as described in 



^The uncertainty due to using a different scheme may be estimated to be of the order of the 0{ag) correction, 
which is about 3 J2 ^qi^sic^s/''^)^- We may take ra = —128 |^ which leads to an uncertainty much smaller than 
that given in Table 1. 

^We thank Thomas Teubner for valuable discussions concerning the perturbative QCD contribution. 
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a reliable prediction for R{s') in the region 1.9 < < 2.8 GeV, then the contribution from 
this interval would be (13.18 ± 0.09) x 10"^. This is essentially the same as the "inclusive" 
contribution, but with a much smaller error. The "exclusive" contribution is a little larger on 
account of the necessity to smoothly match up to the "exclusive" data at ^/s = 1.9 GeV. 

Table 1: Contributions to Aahadl^I) x 10^ of (^) coming from the different ^/s' intervals. 
The alternative values in the round brackets use the summation of exclusive channels as the 
contribution from the region 1.46-1.9 GeV, rather than that from the inclusive data for R{s'). 
Perturbative QCD is used to evaluate the contributions in the intervals 2.8 < y/s' < 3.74 and 
y/s' > 5 GeV. The errors on these contributions are described in the text. 



interval 
(GeV) 


Contribution to 
Atthad X 10^ 


2m^ - 1.46" 
1.46 - 1.9 

1.9 - 2.8 

2.8 - 3.74 
3.74 - 5 

5 - oo 

u, 0, T/^'s, T's 


± { 0.60' } 

8.66 ± 0.60^= 
(10.32 ± 1.06'') 

13.24 ± 0.88^= 
(13.88 ±0.88) 

9.73 ± 0.05'^ 

15.02 ± 0.49 
169.97 ± 0.64^^ 

18.79 ± 0.58 


A«[>a^d X 10^ 


273.82 ± 1.97 
(276.12 ±2.20) 




128.978 ±0.027 
(128.946 ±0.030) 



" The upper (lower) error corresponds to the 2tt (remaining) exclusive channels. 

b,c,d Effors with identical superscripts are added linearly. The remaining errors are added in quadrature. 

For the pQCD contribution we take the mass of the charm quark = 1.46 GeV, and the scale fj? of the QCD coupling 05 to 
be s'. 

To conclude, we see that the introduction of new data has improved the error of Aahad? and 
moved the "inclusive" and "exclusive" determinations a little closer together (compare (^ and 
d^) with the values of Aohad in Table 1). Second, the new BES-II data appear to join more 
smoothly to the inclusive data than the exclusive measurements. However the selection is not 
conclusive and there remains a major residual uncertainty in the 1.5 < ^/s < 1.9 GeV interval, 
which is reflected by the two different results for Acthad listed in Table 1. However if we take 
the favoured "inclusive" result 

Aa[f2j = (273.82 ± 1.97) x 10'^ (10) 

then 

a-\Ml) = 128.978 ±0.027. (11) 
Using this value, the latest profile, obtained^ by the LEP and SLD Electroweak Working 
^We thank Martin Griinewald for making this plot. 
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Group p, for different values of the mass of the (Standard Model) Higgs boson, is shown 
by the dashed line in Fig. 2. We see that it produces a significant increase in the predicted 
mass of the Higgs, with the preferred value lying close to the LEP2 bound. The new 
profile accommodates the LEP2 bound on the mass more comfortably. Note that even the 
determination using the exclusive data up to = 1.9 GeV, now predicts a lower value of 

Aa[f2j = (276.12 ±2.20) x 10-^ (12) 

and a higher Higgs mass than before, with the corresponding profile curve being nearer to 
the dashed curve than the original continuous curve in Fig. 2. 

We emphasize that these standard profiles versus Mu do not include the uncertainty due 
to Acthad- Clearly when confidence limits are placed on the Higgs mass, it is important that 
the uncertainty due to Aohad is included. 

After the completion of this work we were made aware of a recent very preliminary deter- 
mination p5[] , 

Aa£Jj = (275.5 ±4.6) x 10-^ (13) 

which also incorporates the BES-II data. The details of the calculation are not available, but 
from the choice of input for R in the region i/s < 2 GeV we would anticipate a value of Aohad 
somewhat closer to our "exclusive" determination, given in ([12|) . Therefore, if we assume the 
contributions of the other regions are the same, the two results (|T^) and (|13|) are in good 
agreement with each other. 
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Figure 1: A plot of the ratio R{s) versus ^/s in the most sensitive ^/s interval. Up to ^/s = 
2.125 GeV we show by continuous lines the upper and lower bounds of the sum of the exclusive 
channels. However in the "exclusive" analysis we only use these data up to y/s = 1.9 GeV. 
The value quoted in Table 1 corresponds to the average of the bounds. Above y/s = 1.46 GeV 
we show the inclusive measurements of R{s), together with the interpolation (dashed curve) 
used in (H). For 2.8 < < 3.74 perturbative QCD is used to evaluate R; the central QCD 
prediction is shown by the continuous curve. In the region 1.46 < ^/s < 1.9 we used, in turn, 
the inclusive and exclusive data to evaluate Aohad, see Table 1. 
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Figure 2: The latest profile versus the Standard Model Higgs mass obtained by the LEP and 
SLD Electroweak Working Group |2^, which includes the latest preliminary experimental 
data, using Aa^ad of (H) compared to that obtained using (|T0|) , shown by continuous and dashed 
curves respectively. Note that the profiles do not include the uncertainty of Acthad- The 
shaded region to the left is excluded by searches for the Higgs boson at LEP2. 
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